Entanglement Amplification in the Non-Perturbative Dynamics of Modular Quantum Systems 
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We analyze the conditions for entanglement amplification between distant and not directly interacting quan- 
tum objects by their common coupling to media with static modular structure and subject to a local (single-bond) 
quenched dynamics. We show that in the non-perturbative regime of the dynamics the initial end-to-end entan- 
glement is strongly amplified and, moreover, can be distributed efficiently between distant objects. Due to its 
intrinsic local and non-perturbative nature the dynamics is fast and robust against thermal fluctuations, and its 
control is undemanding. We show that the origin of entanglement amplification lies in the interference of the 
ground state and at most one of the low-lying energy eigenstates. The scheme can be generalized to provide a 
fast and efficient router for generating entanglement between simultaneous multiple users. 

PACS numbers: 03.67.-a, 03.67.Hk, 03.65.Ud, 75.10.Pq, 75.20.Hr 



Realizing a sizable and stable entanglement between distant 
qubits is a crucial stage in performing fundamental tasks of 
quantum information and communication . In the ground 
state of noncritical strongly correlated systems with short- 
range interactions the 2-points correlation functions and thus 
the entanglement between spins can be proven to vanish expo- 
nentially with the distance (2]. To overcome this problem, a 
possibility is to mediate indirect interactions between distant 
qubits by a suitable quantum many-body medium. For in- 
stance, one can exploit impurities weakly coupled to the ends 
of a spin chain that mediates an indirect interaction and cre- 
ates a strong long-distance entanglement between them fl. 
Indeed, it has been shown that a relatively large family of 
strongly correlated many-body systems allows for this pos- 
sibility J3-0- However, for most of these systems the energy 
gap is exponentially decreasing in the size of the system, so 
that even for short chains the mechanism becomes thermally 
unstable: very small thermal fluctuations are sufficient to mix 
the ground state with higher energy eigenstates and suppress 
the entanglement between the end impurities. One can intro- 
duce systems with interaction patterns such that the gap closes 
algebraically with the distance [4]; in this case the thermal 
stability improves at the price of letting the end-to-end entan- 
glement become weakly decreasing with the size of the sys- 
tem flllt]. Even if all these results can be extended to higher 
dimensional systems |01 this appears to be the unavoidable 
limit of a a purely static, ground-state approach JH 

On the other hand, it is well known that properly tailored 
time evolutions can propagate entanglement through many- 
body systems @1. In particular, global ifioll or local quan- 
tum quenches 11 ll4l4ll can create long-distance entanglement. 



Apart from certain perturbative proposals that suffer from 



a very slow convergence [11 11—11311 - the dynamical schemes 



in general do not require the weak end coupling assump- 



tion 
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15H and hence thermal instability is not as 
dramatic issue as in the case of static schemes. Rather, the 
price to be payed is the need for a very accurate tuning of the 
times at which optimal entanglement is generated and thus a 
much more elaborate control on the system. In conclusion, 



entangling distant qubits can be realized either via static or 
dynamic approaches, however, either at the price of a strong 
thermal instability or an excessively demanding control. 

In the present work we introduce a scheme that combines 
the modular-static and the quench-dynamic approaches to en- 
tanglement generation and distribution among distant qubits. 
The proposed mixed scheme strongly reduces the drawbacks 
of both approaches and realizes a novel mechanism of en- 
tanglement amplification starting from weakly entangled in- 
puts. In our model, a quantum spin chain is split into two 
elementary bulk modules. Two end impurities are attached to 
each module with couplings of arbitrary strength. This initial 
static configuration is then evolved through a sudden quench 
of the bond connecting the two modules. We show that the 
generated end-to-end entanglement across the entire system 
is always larger than the initial entanglement at the ends of 
each module and that this entanglement amplification can al- 
ways be achieved for all parameter ranges. As the length of 
each module is half the total size of the system and the im- 
purity couplings are non perturbative, the energy gap remains 
sizable even for very long chains and stability against ther- 
mal fluctuations is assured even at moderately high tempera- 
tures. Moreover, by exploiting just a single bond quench for 
inducing a nontrivial dynamics, the required control is mini- 
mal. This mechanism can demonstrate an entanglement router 
which unlike previous proposals neither needs AC-control of 
the couplings llql nor the presence of both ferromagnetic 
and anti-ferromagnetic couplings simultaneously [|17]. Fur- 
thermore, the entanglement amplification mechanism finds a 
clear physical explanation in the quantum interference be- 
tween only two eigenvectors (ground and one of the low-lying 
excited energy eigenstates) of the final Hamiltonian. 

The model. We consider two independent one-dimensional 
arrays of qubits, the modules. Each module is constituted by 
a bulk of spins coupled at both ends with two impurities. The 
interaction Hamiltonian for each module reads: 
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ule has even length Nk and J' k < J so that the spins in the 
bulk tend to entangle between themselves and consequently 
the two end impurities are forced to get entangled From 
the ground state of the XXZ Hamiltonian Hk, one can com- 
pute the initial reduced density matrix of the two impurities 
by tracing out the spins of the bulk. Due to the symmetries of 
the system, the reduced state is diagonal in the Bell basis 
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FIG. 1 : (Color online) (a) A schematic picture of two independent 
modules and their couplings to the respective end impurities, (b) On- 
set of a quenched interaction dynamics between the two modules by 
switching on instantaneously a strong bond Ji between the impuri- 
ties at the module-module boundary. 
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FIG. 2: (Color online) Energy gap of a single XX module as a 
function of: (a) the coupling J' to the end impurities for different 
sizes Nk of the module; (b) the size Nk for different values of J'. 
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where k =L, R denotes the left or right module, respectively 
with N k = N U N R spins, {x[ k \Y^ k \z[ k) } are the Pauli spin 
operators at site i in the fc-th module, J > is the exchange 
spin-spin coupling strength, J' k > specifies the coupling to 
the end impurities, and 6 is the interaction anisotropy along 
the z direction. A schematic picture of this system is shown in 
Fig- [11(a) (note the mirror inversion of spin numbering in each 
module). The spin- 1/2 XXZ Hamiltonian Eq. (Q~|) for each 
module has a rich zero-temperature quantum phase diagram; 
in particular, for — 1 < S < 1, the system is in the so-called 
gapless XY antiferromagnetic phase that admits a nondegen- 
erate, highly entangled ground state. In the two fully isotropic 
limits one recovers the relevant cases of the XX Hamiltonian 
(6 = 0) and the Heisenberg (XXX) Hamiltonian (6=1). 

Static end-to-end entanglement in modules. Let us first dis- 
cuss the entanglement between the two end impurities in the 
ground state of a single module. The model in Eq. <£XJ exhibits 
long-distance entanglement in the ground state when the mod- 
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where is the singlet state, P k is the singlet fraction, \B a ) 
(for a = x,y, z) are the other Bell states which can be ob- 
tained by applying the a Pauli operator on one part of the 
singlet, and the P£"s account for their contributions. For a 
fixed length Nk, one can increase P% arbitrarily and create 
large entanglement between the two impurities by decreasing 
JL However, by decreasing JL the energy gap of the system 



also decreases, at best algebraically |4|,|5|]; when the thermal 
energy ksT becomes comparable with the gap, the state of 
the system becomes a mixed thermal state with no entangle- 
ment between the end impurities. In Fig. [2 a) the energy gap 
A is plotted as a function of J' k for a XX module of length 
Nk which clearly shows an exponential decay for small J' k 's. 
Of course, the gap also decays with the size Nk as shown in 
Fig. Eft). 

Dynamical entanglement generation in a modular quantum 
sjrin chain. Given the initial static situation with the two mod- 
ules in their respective ground states, we introduce a quench 
dynamics between them. The two even-sized modules of 
lengths Nl and Nr are described via the Hamiltonians Hl 
and Hr introduced in Eq. (|T). The exchange coupling J and 
anisotropic parameter 6 are assumed to be identical in both 
modules, while the couplings with the end impurities are re- 
spectively J' L and J' R . The initial ground state of a chain of 
total lengt A^ = Nl + Nr formed by the two noninteracting 
modules is obviously the tensor product of the ground states 
of the two subsystems: \iP(0)) = \GSl)®\GSr). Ati = 0one 
switches on the interaction between the two modules: 
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where Jj is the bond that couples the two modules. The total 
Hamiltonian of the system at t > becomes, as schematicly 
represented in Fig.[TJb), Ht = Hl+Hr+Hj. The initial prod- 
uct state is not an eigenvector of Ht and the system evolves 



according to \tp(t)) ■ 
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two end impurities at time t is 
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We use relative entropy of entanglement [18] as an operational 
measure to quantify the entanglement content of the two end- 
ing impurities at sites 1l and 1r defined by 



E lL , lR {t)=l-H{P s 0Ut {t)) 
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where H(x) = — xlog 9 x— (1 — x) log 9 (l —x) . This is of course re- 
lated to the concurrence IU9I1 of the two impurities Cx Lt \ R {t) 
via E lL . lR (t) = l-H((C lLtlR (t) + l)/2) and always give a 
lower value with respect to concurrence. 

Perturbative regime. We first study the case of J' L and J' R 
both sufficiently small so that the pairs of end impurities in 
each module are initially highly entangled (strong initial end- 
to-end entanglement in both modules). In this situation, in 
each module the impurities are effectively decoupled from the 
rest of the system and resorting to the Schrieffer- Wolff trans- 
formation Il20ll one obtains an effective interaction Hamilto- 
nian for the impurities in each module 

H eff = fitf&A + Y l Y k + 5Z t Z k) . (5) 

where the effective coupling is linear in the energy gap: 
J^ff = Afe/4 with k = L,R. Fixing the inter-module inter- 
action bond at Jj = + Jfft, as shown in Fig. 02b), one 
has that E 1l . 1r (t) = l-H ( 5 ~ 3c ° 8 s(4Jjt) ) . Therefore, at the 

optimal time t op t = jjj, the maximal long-distance entan- 
glement is established between the ending sites. Compared 
to the static case with a single long module of length N with 
the same entanglement between the end impurities 1 l and 1^, 
the dynamical mechanism in the perturbative regime assures 
a larger thermal stability "per se". If the entangling time t op t 
is engineered so to be much smaller than the thermalization 
time, then the thermal effects are fully determined by the ther- 
mal initial state associated to the energy gap of each module, 
which is always above and can be made much larger than the 
gap of the two combined in a single one of size N = Nl+Nr. 
In the dynamic case one can exploit larger impurity couplings 
J' k for each module in comparison to the static case with a sin- 
gle module of size N, thus increasing the gap, as illustrated in 
Fig. 12 a). In this way thermal instability is ameliorated but 
not fixed; moreover, the perturbative nature of the coupling 
J 1 implies a very slow dynamics. The power of the scheme 
shows only in the non-perturbative regime. 

Non-perturbative regime: Entanglement amplification. 
When the strength of the coupling to the impurities becomes 
comparable to the interaction energies in the bulk, the ini- 
tial end-to-end entanglement for each module is strongly sup- 
pressed, the reduction to the effective Hamiltonians Eq. @ 
is no longer justified, and Jj cannot be determined analyti- 
cally as in the perturbative regime. To proceed, one can ob- 
serve that in practice the affordable time for the entangling 
dynamics is ultimately bounded by the decoherence rates. In- 
troducing such upper bounds for the optimal entangling time, 
we define an optimization problem for the largest attainable 
end-to-end entanglement that depends on two free parame- 
ters, i.e. the static coupling to the end impurities J' = J' L = J' R 
and the quench-dynamic bond J/ between the two modules. 
In Fig. [3 a) the end-to-end entanglement Ex L ^ R (t) is plotted 
as a function of time both for the Heisenberg and the XX 
Hamiltonians for different values of J 1 and Jj. As reported 
in Fig. |3ja), the time evolution after the quench generates 
a large end-to-end entanglement, strongly amplified respect 




FIG. 3: (Color online) (a) End-to-end entanglement as a function of 
time for Nl=Nr=W, J'— 0.5 and J/=0.65. (b) E max as a function 
of N. (c) Optimal time t op ± versus length N. (d) Optimal impurity 
coupling J' opt for XXX (open blue circles) and XX (green open di- 
amonds) and optimal interaction coupling Jj pt for XXX (filled red 
circles) and XX (filled brown diamonds). 



to the initial one in each module. The optimal entangling 
time t op t is then the earliest time at which the end-to-end 
entanglement peaks, defining its maximum attainable value: 
E max = E(t opt ). In Fig. Ob) we plot E max as a function 
of the total chain's size N when both J' and Jj are tuned to 
their optimal values. One has that for long enough chains both 
the XX and Heisenberg Hamiltonians perform equally well. 
However, if we look at Fig. |3jc), where t opt is reported as a 
function of N, one can see that the Heisenberg chain gener- 
ates a faster dynamics with earlier peaks, an important added 
value in order to minimize the effects of decoherence. Fi- 
nally, in Fig.[3jd) we report the optimal values J' opt and J° pt , 
respectively of the coupling to the end impurities and the in- 
teraction coupling in the bulk, as functions of N. Remark- 
ably, as Fig. 0d) clearly shows, one finds that the optimal 
couplings decrease very slowly by increasing the size N, im- 
plying the onset of a fast dynamics and the permanence of 
the non-perturbative regime even for very long chains. Fur- 
thermore, larger couplings in Heisenberg chains in compare 
to XX ones results in a faster dynamics and provides higher 
energy gap and thus higher thermal stability. 

To show how entanglement is actually amplified by the dy- 
namical process we compare the initial static end-to-end en- 
tanglement Ei L N t (0) in the ground state of a single mod- 
ule of length Nl = N/2 and the output dynamical end-to- 
end entanglement Ei L ^ R (t opt ) at the optimal entangling time 
across two interacting modules forming a chain of total length 
Nl+Nr = N. Fig.|4{a) shows that the initial ground-state end- 
to-end entanglement is always amplified for the entire range 
of couplings J' to the end impurities. A remarkable rebound 
of Ei Lt i B (t op t) occurs at a point of non-analyticity that sep- 
arates the perturbative from the non-perturbative regime, cor- 
responding to the onset of enhanced dynamical amplification 
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FIG. 4: (Color online) (a) The initial, static end-to-end entangle- 
ment Ei L ,isr L (0) (or Ei r ,n r (0)) in each single module of length 
Nl = Nr — 10 together with the maximal end-to-end dynamical 
entanglement Ei L ,i R (t op t for the entire two-module chain of length 
N = Nl + Nr as functions of the coupling J' for 5=1. (b) The 
optimal time t opt at which entanglement peaks, as a function of J'. 



against pure decay. In order to characterize and understand 
these two different regimes and the transition between them 
we have studied the entangling time t opt as a function of the 
coupling J' to the end impurities. Fig.|4](b) illustrates that in 
the perturbative regime t opt increases exponentially with de- 
creasing J', while in the non-perturbative regime it remains 
essentially fiat at very small values, guaranteeing, for appro- 
priately chosen values of J', that maximal amplification oc- 
curs well in advance of the effects of decoherence. 

Origin of amplification: Excitation spectrum and quantum 
interference. To understand the physical mechanism under- 
lying entanglement amplification in the non-perturbative dy- 
namics of modular many -body systems we must take notice of 

the fact that|V(t))=E„=i c « e " " where, is 
the n-th eigenvalue (eigenvector) of Ht and c n = (-B n |^(0)). 
When the values of the Hamiltonian parameters J' and Jj are 
far from their optimal values, many of c„'s are significantly 
different from zero. At values close to the optimal set the situ- 
ation becomes radically different, as reported in Figs. |5J a ) an d 
(b), where the squared amplitudes \c n \ 2 are plotted as func- 
tions of the energy difference between the corresponding n-th 
eigenstate and the ground state of Ht, respectively for the 
XX (5 = 0, Fig. [2a)) and the Heisenberg Hamiltonian (6=1, 
Fig. |3Ja)). In both cases the initial state is essentially pro- 
jected onto only two eigenstates of Ht, the ground state and 
just one of the first few low-lying excited states. Therefore the 
evolution of the initial state at time t under the action of Ht is 
essentially due to the relative phase factor <f>(t) = exp(— iwgt), 
where ojs is the the difference of the eigenvalues of the two 
eigenstates as indicated in Figs. [5ja) and (b). As shown in 
Figs. |3c) and(d), the maximum of the entanglement ampli- 
fication is reached when <f)(t) ~ — 1. Therefore the dy- 
namic amplification of long-distance entanglement is essen- 
tially due to a constructive interference between the two eigen- 
state while the fact that the maximum is reached around and 
not exactly at (j>(t) = — 1 accounts for the projections on the 
remaining part of the spectrum: although strongly suppressed, 
they are not exactly vanishing. As we move to values of the 
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FIG. 5: (Color online) Upper panels: |c^| versus AE/J = (E n — 
Eo) /J for all eigenvectors of Ht for two chains with Nl = Nr = 4, 
J' L = J'r = 0.5J, Ji — 0.75J where J is the bulk amplitude for (a) 
5 = 0; (b) 5 = 1. Lower panels: Ei Ll i R vs. time for the two chains 
with parameters given above and (c) 5 = 0; (d) 5 = 1. The red 
dashed lines stand for the initial end to end entanglement in each 
single chain. The vertical gray lines signal the times at which the 
phase (/>(£) =exp(— iugt) =— 1. 

Hamiltonian parameters away from the optimal ones and/or 
the dimension of the two modules is increased, the relative 
weight of the other eigenstates becomes more relevant reduc- 
ing the maximum value reached by Ei L! \ R . 

Conclusions & Outlook: Entanglement router. By consid- 
ering many-body media with static structures of bulk mod- 
ules and by engineering suitably quenched local interaction 
dynamics between different modules, we have introduced 
a method for the generation and amplification of entangle- 
ment between distant and non-interacting quantum objects. 
The method is intrinsically non-perturbative and does not re- 
quire sophisticated controls of the system dynamics. Indeed, 
we showed that the end-to-end entanglement initially present 
in the ground state of a modular many-body system is am- 
plified by the quenched dynamics whenever all the values 
of the Hamiltonian parameters are far from the perturbative 
regime. The ensuing entangling dynamics is fast and robust 
against thermal fluctuations. The occurrence of such dynam- 
ical mechanism of entanglement amplification finds a sim- 
ple and beautiful explanation in the constructive interference 
of only two energy eigenstates of the driving Hamiltonian. 
This opens potentially many new perspective for studying en- 
tanglement generation, distribution, and manipulation across 
arbitrary distances in a many-body systems. Our minimal 
control strategy applies immediately to the demonstration of 
an entanglement router capable of distributing entanglement 
simultaneously between multiple users. Indeed, consider a 
setup in which every user controls one impurity in a mod- 
ule which extends from its position to a common dispatching 
center where the other impurities are controlled. In the dis- 
patching center, one can switch on the interaction between 
any pair of impurities of different modules and switch on the 
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quenched entangling dynamics in those particular chains, thus 
creating entanglement directly between the two users. The 
main advantage of such an entanglement router is the mini- 
mal control needed (a single bond quench) that, at variance 
with previous router proposals, neither requires AC-control of 
the couplings fl6ll nor the introduction of many different cou- 
plings El. 
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